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ABSTRACT 

A great deal of interest has been recently raised in connection with the possibility that soft 7- 
ray repeaters (SGRs) and anomalous X-ray pulsars (AXPs) contain magnetars, young neutron 
stars endowed with mag netic fields ^ 10 14 G. In this paper we calculate thermal spectra from 
ultra-magnetized neutron stars for values of the luminosity and magnetic field believed to be 
relevant to SGRs and AXPs. Emergent spectra are found to be very close to a blackbody at the 
star effective temperature and exhibit a distinctive absorption feature at the proton cyclotron 
energy E c>p ~ 0.63(5/10 14 G) keV. The proton cyclotron features (PCFs) are conspicuous 
(equivalent width of up to many hundreds eV) and relatively broad (AE/E ~ 0.05 — 0.2). The 
detection of the PCFs is well within the capabilities of present X-ray spectrometers, like the 
HETGS and METGS on board Chandra. Their observation might provide decisive evidence in 
favor of the existence of magnetars. 

Subject headings: radiative transfer — stars: neutron — X-ray: stars 



1. INTRODUCTION 

Over the last few years increasing observational ev- 
idence has gathered in favor of the existence of ultra- 
magnetized neutron stars (NSs) with surface field 
B ;> 10 14 G. The existence of these magnetars was 
first proposed nearly 10 years ago by Thompson & 
Duncan ( 1993| ) who showed that, during the core col- 
lapse, convective motions can strongly amplify the 
seed magnetic field. 

Magnetars are unique among compact stars, in that 
their luminosity should be provided by the dissipa- 
tion of magnetic energy. Like conventional radio- 
pulsars, however, they will also lose rotational energy 
according to the magnetic dipole formula. Owing to 
their high magnetic field, magnetars are expected to 
spin down at a high rate, P m 10- n (5/10 14 G) 2 /P 
ss -1 . It was the detection of pulsations with a sec- 
ular spin-down in this range in two soft 7-repeaters 
(SGRs; Kouveliotou et al. |1998| ; |199S| ) that for the 
first time suggested the association of these sources 
with ultra-magnetized NSs. Besides their bursting 
activity, SGRs show also persistent X-ray emission 
with luminosities L ~ 10 34 — 10 36 ergs -1 . X-ray 
spectra are usually power-law like, but in one case 
(SGR 1900+14) clear evidence was found for a ther- 



mal spectral component at kT ~ 0.5 keV (Woods 
et al. 1999] ) . This is believed to originate from the 
NS surface radiating away the dissipated magnetic 
energy (see e.g. Thompson 2000 for a recent review). 

Magnetars have been also invoked to explain the 
characteristics of another enigmatic class of Galac- 
tic sources, the so-called anomalous X-ray pulsars 
(AXPs). These pulsars share properties that are 
much at variance with those of canonical X-ray pul- 
sators (Mereghetti & Stella 1995j ). Among these 
are: i) a narrow interval of pulse periods (6-12 s); 
ii) very soft X-ray spectra, often with a blackbody- 
like component at kT ~ 0.4-0.6 keV; iii) relatively 
low X-ray luminosities (~ 10 35 ergs -1 ); iv) spin- 
down trends at relatively large and stable rates (~ 
0.6-30xl0 -12 ss -1 ); v) absence of massive compan- 
ion stars and vi) association with supernova remnants 
(see e.g. Mereghetti 2000 for a review). 

Though it is now widely recognized that AXPs are 
rotating NSs, two different mechanisms for powering 
their X-ray emission have been discussed: accretion 
(from a very low mass companion, fall back or a fossil 
disk, Van Paradijs, Taam, & Van den Heuvel |1995|; 



Chatterjee, Hernquist, & Narayan |2000| ) and mag- 
netic energy dissipation in ultra- magnetized, isolated 



1 Mullard Space Science Laboratory, University College London, Holmbury St. Mary, Dorking, Surrey, RH5 6NT, UK; e-mail: 
sz@mssl.ucl.ac.uk 

2 Dipartimento di Fisica, Universita di Padova, Via Marzolo 8, 1-35131 Padova, Italy; e-mail: turolla@pd.infn.it 
3 Osservatorio Astronomico di Roma, via dell'Osservatorio 2, 1-00040 Monte Porzio Catone, Roma, Italy; e-mail: 
stella@coma.mporzio.astro.it 

4 Dipartimento di Scienze, Universita delHnsubria, Via Valleggio 11, Como, Italy; e-mail: treves@mib.infn.it 



2 



PROTON CYCLOTRON FEATURES 



NSs. If AXPs are young, hot NSs giving off thermal 
radiation, magnetic field decay should provide crustal 
heating for a long enough time. Possible decay sce- 
narios have been analyzed by Colpi, Geppert & Page 
( pOOOj ) and led to the conclusion that AXPs are ~ 10 4 
yr old with a present 5-field of ~ 10 14 G. The latter 
figure is in the range of that inferred from SGRs spin- 
down rates. The many similarities between AXPs 
and SGRs strengthen the idea that the two classes 
of sources are powered by the same mechanism, i.e. 



magnetic field dissipation (e.g. Thompson 2000). 

The transfer of radiation in strongly magnetized 
atmospheres has been widely investigated in connec- 
tion with X-ray emission from cooling and accreting 
neutron stars (see e.g. Shibanov et al. 1992; Ra- 
jagopal, & Romani 1996; Pavlov et al. 11996 ; Zane 



Turolla & Treves |200C , herafter ZTT). Up to now, 
however, little attention was given to emission from 
ultra-magnetized neutron stars with surface fields 
J> 10 13 G. Detailed transfer calculations to assess the 
emergent spectrum from atmospheres around ultra- 
magnetized NSs are therefore needed. Model fits 
to the X-ray spectra of AXPs and SGRs indicate 
the presence of a power-law, high-energy component. 
In addition to this, thermal emission has been de- 
tected in all A XPs b ut one, and in SGR 1900+14 
(see Mereghetti 2000| and references therein). 

In this paper we present numerical models of ther- 
mal emission from NSs for a range of luminosity, 
L ~ 10 34 — 10 36 ergs -1 , and magnetic field, B ~ 
10 14 - 10 15 G, which are thought typical of SGRs 
and AXPs. Emergent spectra, computed solving the 
full (angle and energy dependent) transfer equations, 
are nearly Planckian, have negligible hardening, and 
exhibit a strong absorption feature at the proton 
cyclotron resonance, E CtP ~ 0.63i/g(B /W 14 G) keV, 
where yc = \A — 2GM/c 2 R is the gravitational red- 
shift factor (M and R are the star mass and radius, 
respectively). The predicted equivalent width of the 
proton cyclotron features (PCFs) is in the range 0.03- 
1 keV and the relative width FWZI/E CjP « 0.7 (here 
FWZI is the full width at zero intensity). The detec- 
tion of PCFs is definitely within the range of Chandra 
spectrometers. Their actual observation in the X-ray 
spectra of SGRs and/or AXPs would provide a di- 
rect measurement of the surface magnetic field and a 
remarkable confirmation of the magnetar model. 

2. MODEL ATMOSPHERES 

In this paper we follow closely the approach pre- 
sented in ZTT, to which we refer for details. In partic- 
ular we consider a magnetized, non-degenerate, pure 
hydrogen plasma in thermal equilibrium at tempera- 
ture T in which the dominant radiative processes are 
free- free emission/absorption and Thomson scatter- 



ing; non-conservative scattering has negligible effect 
on the transfer of radiation but becomes important 
in the energy balance of the outer layers (see below). 
We approximate the atmosphere as a plane-parallel 
slab with the magnetic field B parallel to the nor- 
mal n (the x-axis). This is reasonable for magne- 
tized accretors, but it becomes questionable in treat- 
ing emission from the entire star surface, as in the 
case of cooling neutron stars. The effects induced on 
the emerging spectrum by including surface elements 
away from the magnetic poles are therefore discussed 
in section ||. Although we are mainly interested in 
thermal emission from magnetars, ultra-magnetized 
NSs might also accrete the interstellar medium, as 
recently suggested by Rutledge ( [200 1| ). In the light 
of this, and for the sake of completeness, we consider 
here both the case of accreting and cooling atmo- 
spheres. 

For photon frequencies in the range u> p < uj < ui c ,e, 
where u> p and uj Cte = eB/m e c are the plasma and the 
electron cyclotron frequency respectively, the semi- 
transverse approximation can be used and the trans- 
fer of radiation through the magneto-active medium 
is described in terms of the two normal modes (ZTT 
and references therein). Atmospheric models are con- 
structed by solving the transfer equations for the nor- 
mal modes coupled to the hydrostatic equilibrium 
and the energy equation. Since the luminosity L is 
well below the Eddington limit, the hydrostatic bal- 
ance gives immediately the density as a function of 
the scattering depth r = n es pdz 



GMm r 



P 



2y 2 G R 2 K es fcT(r) 



(1) 



(here K es = 0.4 cm 2 g _1 is the electron scattering 
opacity for an unmagnetized medium; see ZTT). The 
energy balance simply states that the net radiative 
cooling must equate the energy Wjj released in the 
atmosphere by non-radiative processes, like the bom- 
bardment of impinging ions; for cooling atmospheres 
Wh = 0. Including Compton heating/cooling the 
energy equation becomes 



aT 4 



5>«£^ } + (r-A) 



c 



CK P 



where U^ ' is the radiation energy density of mode i, 
kp, ham are defined in strict analogy with the Planck 
and absorption mean opacities in the unmagnetized 
case and (r — A)c is the Compton heating/cooling 
rate in a magnetized plasma (see again ZTT). 
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3. RESULTS 

Following the method outlined in the previous sec- 
tion, several model atmospheres have been computed 
for different values of the magnetic field strength 
and luminosity. As expected (see ZTT; Treves et 
al. |2000| ), cooling and accreting atmospheres turn 
out to produce essentially the same spectrum at X- 
ray energies, so we do not make further distinction 
between them. It should be kept in mind, however, 
that a cooling model characterized by a luminosity 
L is equivalent to an accreting one with luminosity 
JaL < L, where /a is the fraction of the emitting 
area. All models have been calculated for M = 1M & , 
R = 9 km and are labeled by the value of the lu- 
minosity measured by an observer at infinity, assum- 
ing emission from the entire surface (i.e. Ja = 1)- 
We have explored the ranges 10 13 G < B < 10 15 G 
and 10 34 ergs _1 L <^ 10 36 ergs _1 which are be- 
lieved typical of magnetars. The main parameters 
of the computed models are summarized in table |l] 
and emerging spectra are shown in figure |. In all 
cases, with the exception of model A4, the spectra 
are nearly Planckian in shape and show a small hard- 
ening with respect to the blackbody at the star ef- 
fective temperature T e ff. As it was found in previ- 
ous investigations (e.g. Shibanov et al. 1992) ; ZTT), 
the spectral hardening decreases with increasing field 
strength and almost vanishes for B ^ 10 14 G. The 
most prominent spectral signature is the absorption 
feature at the proton cyclotron resonance E cp ~ 
0.63?/g ! (^/10 14 G) keV which falls in the soft-medium 
X-rays for B ~ 10 -10 15 G; the electron cyclotron 
line is at ~ 1.2(5/10 14 G) MeV. The line equivalent 
width, EW, and the full width at zero intensity rel- 
ative to the line center, AE/E = FWZI/E, are re- 
ported in table []]. Two main effects contribute to 
the line profile: the intrinsic resonance in the mag- 
netic absorption coefficients (this essentially produces 
Fraunhofer lines) and mode crossing at the mode 
collapse points. The latter is amplified when col- 
lapse points induced by vacuum effects fall near the 
resonance and in the photospheric region (see also 
Shibanov et al. 1992). The large modification in the 



continuum in model A4, for instance, is partially due 
to mode-crossing which, for these values of the pa- 
rameters, occurs close to the Wien peak. 

As discussed earlier, the assumption of a constant 
magnetic field breaks down when emission comes 
from the entire star surface. Even for a simple dipo- 
lar field, B changes in both magnitude and direction 
along the surface and this will produce a broaden- 
ing of the PCF. Moreover, large -B-fields introduce 
a meridional temperature variation that makes the 
emergent spectrum dependent of the viewing angle. 
A thorough modeling would require a genuine two- 



dimensional transfer calculation, which is still beyond 
the capabilities of present numerical codes. However, 
an estimate of the expected broadening can be ob- 
tained with a simple, approximated computation. To 
do this, we assume a relativistic magnetic dipole in 
the Schwarzschild space-time, aligned with the rota- 
tional axis (see e.g. Pavlov & Zavlin [2000) ) 



B 



B u 



(4 - / 2 ) cos 2 X + f 



1/2 



(3) 



where B p is the field strength at the magnetic pole 
(x = 0) an d / — 1.1 accounts for gravitational ef- 
fects. The surface temperature profile is taken from 
Possenti, Mereghetti & Colpi (|1996[) 



'eff 



K+{A-K) cos 2 x 



1/4 



[1-0.47(1 -A")] (1 + 3 cos 2 X ) 



(4) 



where K is the ratio of the coefficients of thermal 
conductivity orthogonal and parallel to the field. The 
exact value of K is not very important: for B ;> 10 13 
G it is K <; 10~ 4 , and beyond this value the temper- 
ature profile depends weakly on K. The flux emitted 
by the surface element at co-latitude x is evaluated 
assuming thermal diffusion for both modes 



dF u 



2 J>(0 dB u [T( X )} . , 
2vr2^— — sm X dx 



dr 



(5) 



where B v is the blackbody function. The diffusion 
coefficients D^' are obtained by angle-averaging the 
reciprocal of the total opacity over all photon direc- 
tions (see Shibanov et al. 1992; ). They depend on 
Hb = n • B/S; here n = r/r is the "local" normal 
to the surface. The total flux is then computed by 
integrating equation (|5|) over the entire star surface, 
with the aid of expression 



1 



4-/ 2 



4/ 2 



(4-/ 2 )cos 2 X + / 2 



(6) 



which relates ^ to x (see again Pavlov & Zavlin 
2000|) . 

Line parameters derived from the previous calcu- 
lation are reported in table ^ for different values of 
the luminosity and of the polar S-field; diffusion-limit 
spectra, as seen at the top of the NS atmosphere, are 
shown in figures g]| For B v = 10 15 G the PCF 
falls in the Wien tail and it is spread over a few 
keVs. To allow a direct comparison with the case 
of a magnetic field parallel to the surface normal, 
we have repeated the same calculation for a plane- 
parallel slab with constant B = B p , neglecting sur- 
face temperature variations. We stress that in both 
cases no mode crossing effects have been accounted 
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Table 1 
Model Parameters 



Model 


B 
10 13 G 


L 

10 34 erg s _1 


pa 

keV 


EW 
keV 


AE/E 


Al 


1 


2 


0.06 


0.01 


0.67 


A2 


1 


70 


0.06 


0.01 


0.47 


A3 


5 


9 


0.32 


0.05 


0.31 


A4 


10 


3 


0.63 






A5 


50 


170 


3.15 


0.10 


0.10 


A6 


100 


180 


6.3 


0.11 


0.06 



Note. — The line energy is not corrected for the gravita- 
tional red-shift; the corresponding value at Earth is lower 
by a factor yc ~ 0.8. 



for, and, to be consistent with the diffusion approxi- 
mation, spectra were computed at optical depth ^ 1. 
Clearly these approximated results may be taken only 
as indicative of the true shape of the PCF in emer- 
gent spectra; it is however worth noticing that the 
line equivalent widths derived from our simple analy- 
sis for B 1 1 n are in rough agreement with those of the 
complete transfer calculation for the same values of B 
and L. As expected, the PCF turns out to be system- 
atically broader when emission comes from the entire 
star surface, typically by 10-20%. Also, the change of 
the field strength with the latitude produces a shift of 
the line centroid toward lower energies by ~ 20-30%. 
Both these effects are quite independent of the val- 
ues of Bp and L (see table |2|). Complex magnetic field 
configurations with contributions from higher multi- 
poles, such as those envisaged in magnetar models 
(Thompson |2000| ) or in the modelling of isolated neu- 
tron stars (e.g. Geminga, Page et al. 1995|) , perhaps 
represent a situation more realistic than the dipolar 
field we have considered. Due to the steep effective 
temperature gradients expected in this case, the real 
properties of PCFs should be intermediate between 
those discussed above. 

4. DISCUSSION 

In this paper we have modeled thermal emission 
from ultra-magnetized NSs, showing that the emer- 
gent continuum is almost indistinguishable from a 
blackbody at T e ff. We point out the existence of a 
strong absorption feature at the proton cyclotron en- 
ergy; the line equivalent width is between 30 eV and 



1 keV and the line center is located at ~ 0.5-5 keV 
for field strengths ~ 10 14 -10 15 G, which are currently 
believed to be typical of SGRs and AXPs. The radia- 
tive transfer calculations presented here corroborate 
the original suggestion by Thompson ( |2000 ) that an 
emission/absorption feature at the proton cyclotron 
resonance should be present in the thermal spectrum 
of magnetars. 

The natural width and the line-splitting of elec- 
tron cyclotron lines were thoroughly investigated 
by Pavlov et al. (1991) by means of a rigorous 
QED treatment. No similar calculation is available 
for protons, but, in any realistic situation, B -C 
B cr it,p = rripC^/he ~ 1.5 x 10 20 G, so protons are 
never relativistic. In this case, the natural width 
of the line associated with the transition between 
the first and the ground Landau levels is rf ~ 
(4/3)am p c 2 (B/B cr it t p) 2 (Pavlov, private communica- 
tion), where a is the fine structure constant and the 
plus/minus sign refers to spin- up *-* spin-up, spin- 
down *-* spin-down transitions respectively. In addi- 
tion to the first order cyclotron resonance Thomp- 
son ( |2000D also suggested the possibility to detect 
the spin-flip line, despite the fact that no detailed 
calculation has been presented. In fact, contrary to 
electrons, which have a small anomalous magnetic 
moment A^/hq ~ a (here /iq is the Bohr magne- 
ton), A^l/^ln — 2.79 for protons, /x/v = ijn v jm^)\iQ 
being the nuclear magneton. This implies that the 
proton spin-flip transition produces a line at 2.792? C)P 
which is always well-separated from the fundamen- 
tal resonance at E C)P . Here we only note that this 
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Table 2 
Line Broadening 



B 






EWn/EWu 

X-J V V J^J 1 J— J v r 




^c,d/^c,\\ 




13 p 




KB V 




WAT" 

Ke V 






1 


0.1 


0.035 


1.17 


0.046 


0.73 


0.84 


1 


1 


0.035 


1.20 


0.046 


0.73 


0.84 


1 


10 


0.035 


1.22 


0.047 


0.75 


0.80 


10 


0.1 


0.35 


1.14 


0.49 


0.78 


0.75 


10 


1 


0.35 


1.12 


0.47 


0.75 


0.86 


10 


10 


0.34 


1.08 


0.46 


0.73 


0.73 


50 


1 


1.74 


1.06 


2.57 


0.81 


0.53 


50 


10 


1.72 


1.06 


2.72 


0.86 


0.63 


50 


100 


1.75 


1.11 


2.49 


0.79 


0.79 



Note. — D = dipolar field; || = constant B field. E c is the line centroid and AEd 
is the FWZI in the case of dipolar field; energies have not been corrected for the 
gravitational red-shift. 



second line should be much weaker and is probably 
beyond the capabilities of the present spectrometers, 
because the rates for the spin-flip transitions in the 
non-relativistic regime should contain an additional 
factor ~ E^p/rripC 2 <C 1 with respect to transitions 
with conserved spin projection (Pavlov, private co- 
munication; see also Melrose & Zheleznyakov 1981 for 
the case of electrons). On the other hand, the fact 
that the lowest Landau level is always non-degenerate 
for protons might prove decisive in discriminating 
between proton and electron absorption features in 
high-resolution spectral observations. In fact, if ab- 
sorption features observed at 0.06 — 6 keV are inter- 
preted as electron cyclotron resonances, the inferred 
magnetic fields should be as low as 10 9 — 10 11 G. For 
values in this range, the interaction with the radia- 
tion field should remove the degeneracy of the lowest 
Landau level, giving rise to a doublet. In the case 
of electrons, even if the two lines are not individu- 
ally resolved (the required resolving power should be 
~ E ce /AE ~ (0.16a) _1 ~ 860) observations with a 
narrow band Bragg polarimeter should show a sharp 
jump when scanning the feature from the blue toward 



the red wing (see Pavlov et al. 1991 for details). 

Finally, we would like to stress that present results 
have been derived under a set of simplifying assump- 
tions, on both the atmospheric properties and the line 
profiles. Atmospheres around young NSs are likely 
to be rich in metals. Unfortunately the absorption 



coefficients for metals in a strong magnetic field are 
poorly known and our assumption of a pure H com- 
position, which is certainly inadequate, appears at 
present the only viable option. Metal lines may well 
be present in magnetars X-ray spectrum. If originat- 
ing in an accreting gas or a corona they are expected 
to be mainly in emission and narrower, unless the 
circumstellar material where the line originates has a 
velocity large enough to produce a (Doppler-) 
broadening comparable to that of the PCF. This is 
not obvious if the NS atmosphere itself is rich in 
heavy elements since energies and shapes of atomic 
transitions are strongly dependent on the local con- 
ditions (see Rajagopal, &; Romani |1996j ; Zavlin et 
al. 1996 ). Strongly magnetized Fe spectra presented 
by Rajagopal, Romani, & Miller ( |1997| ) show signif- 
icant absorption features in the X-rays. A confusion 
between metal lines and PCFs is therefore possible 
and this issue should be carefully addressed when ob- 
served spectra are analyzed. Moreover, for a different 
chemical composition, lower energy cyclotron lines at 
the various ion resonances should also be present, 
since E c i = (Z/A)E CtP . A firm identification of an 
observed feature in terms of a PCF may follow from 
the detection of higher harmonics (already observed 
in some pulsars for the electron cyclotron resonance), 
despite they are expected to be weaker (but see Co- 



sumano et al. 1995 for a counterexample). In the 
normal mode approximation, higher-order harmonics 
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of the electron cyclotron resonance can be accounted 
for in the Gaunt factors (see Pavlov & Panov 1976| ; 
these expressions have been used here and in ZTT), 
but we are unaware of a similar treatment of the ef- 
fective collision frequencies for the higher harmonics 
of the proton resonance. 

Partial H ionization should also be taken into 
account. Although no detailed investigation on 
the ionization properties of hydrogen in ultra-strong 
fields have been carried out, the results of Potekhin, 
Chabrier & Shibanov dl999|) for B £ 10 13 G clearly 
indicate that the neutral fraction increases with the 
field strength. For values of density and temperature 
typical of the inner layers of our atmospheric models, 
the neutral fraction is ~ 0.01 at B = 10 13 G and it 
is bounded to increase at larger B. This will reduce 
the number of free protons, thus making resonant ab- 
sorption less effective. 

Line-broadening mechanisms have not been ac- 
counted for in our present investigation, apart from 
the smearing of the line introduced by the dipo- 
lar variation of the -B-field (see §||). The more im- 
portant broadening could be due to the thermal 
Doppler effect, which produces a width AEth ~ 



Up to now no positive detection of any spectral fea- 
ture, both in emission and in absorption, have been 



.2x1/2 



10 E cp under typical condi- 



E CtP (kT/m p c 

tions. As expected, Doppler-broadening is less im- 
portant by far for protons with respect to electrons 
and AEth "C AE in our models (see table 0). 



reported in SGRs and AXPs (see Mereghetti |200C 
and references therein). The only possible excep- 
tion is the anomalous X-ray pulsar IE 2259+586, for 
which the detection of a cyclotron feature around 5 
keV has been claimed (Iwasawa, Koyama & Halpern 
1992|) but never confirmed by subsequent observa- 
tions (e.g. Corbett et al. 1995). The discrepancy 
between the spin-down age and the age of the SNR 
in this source hints toward a higher value of the field 
in the past, but the large line energy would imply 
a present field ~ 10 15 G, which seems difficult to 
reconcile with the pulsar evolutionary track (Colpi, 
Geppert & Page |2000| ). However, present data do 
not exclude the possibility that in SGRs and AXPs 
the magnetic field is a few 10 G. If this is the case, 
a cyclotron feature around ~ 1 keV may well have 
escaped detection because these sources have never 
been observed with sufficient spectral resolution at 
low energies. 



We are grateful to G.G. Pavlov for several enlight- 
ening comments on the radiative width of the proton 
cyclotron line and the relative importance of the pro- 
ton spin-flip transition. Work partially supported by 
Italian MURST under grant COFIN-98-02154100. 
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Fig. 1.— 
temperature 
evaluated at 



Emergent spectra (solid lines) for the models presented in table 1, together with the blackbody at the NS effective 
(dash-dotted lines). Energies have been red-shifted at Earth while the monochromatic flux (on the vertical axis) is 
the star surface. 






Fig. 2. — Broadening of the proton cyclotron line for B = 10 14 G and different luminosities. From left to right: L = 
10 33 , 10 34 , 10 35 ergs" 1 ; solid lines correspond to a dipolar field, and dash-dotted lines to B normal to the surface. Flux is in 
arbitrary units and energies have not been red-shifted at Earth. 
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Fig. 3.— Same as in figure | for B = 5 x 10 14 G and L = 10 34 , 10 35 , 10 36 ergs" 1 . 
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Fig. 4.— Same as in figure § for B = 10 15 G. 
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